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We make a critical study of two highly-constrained models of supersymmetry — the constrained min- 
imal supersymmetric standard model (cMSSM), and the non- universal Higgs mass model (NUHM) 
— in the light of the 125-126 GeV Higgs boson, the first observation of Bs — >■ ^i^i at the LHCb, and 
the updated B tv branching ratio at BELLE. It turns out that these models are still allowed by 
the experimental data, even if we demand that there be a light stop with mass less than 1.5 TeV. 
The only significant effects of all these constraints are to push the mass of the light stop above 
~ 500 GeV, and to prefer the universal trilinear coupling Aq to be large and negative. We calculate 
the Higgs boson branching ratios to WW, ZZ, tt and 77 in these models and show that improved 
experimental limits on these could put them to the most stringent experimental tests yet. 

PACS numbers: 12.60.Jv, 13.85.Rm, 14.80. Ly, 14.80.Nb 
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The recent discovery at the CERN LHC, of a bo- 
son with mass around 125-126 GeV and properties closely 
resembhng those of the Standard Model (SM) Higgs bo- 
son 0, 0], matches well with a common prediction of 
all supersymmetric models viz., the existence of a 
light Higgs boson with mass in the 100 — 135 GeV range. 
However, as this could be sheer coincidence, only the dis- 
covery of a superpartner of one of the known SM parti- 
cles would really establish the existence of supersymme- 
try (SUSY). Now, the elevation of the light Higgs boson 
mass to a value above Mz « 91 GeV in supersymmet- 
ric models is known to be primarily due to radiative 
corrections involving the top quark (t) and the light and 
heavy stop states {ti,t2). It follows, therefore, that the 
discovery of a light Higgs boson candidate in the upper 
part of the SUSY-allowed range may be considered as a 
hint of the existence of a rather low-lying stop state {ti) 
— perhaps light enough to be discovered at the LHC. We 
cannot, however, be too sure about this, for SUSY mod- 
els can be practically indistinguishable from the SM if we 
take the so-called decoupling limit, where all the super- 
partners are too heavy to be detected at the LHC (or any 
terrestrial experiment conceivable in the near future). 

Faced with this situation, we approach the problem in 
a pragmatic way, viz. we focus on that part of the large 
SUSY parameter space that could support a light stop 
state discoverable at the LHC. Though this does resemble 
the proverbial search under the lamp post, at the present 
juncture, it seems the most reasonable thing to do. The 
question then arises as to how light the ti state must be. 
An exact answer is, of course, a matter of detail, but one 
can form a crude estimate based on the fact that the pro- 
duction of titi pairs at the LHC is dominated by strong 
interactions and hence depends essentially on the mass 
of the ti alone. The variation of the <iii pair-production 
cross-section with m^^ at the 14- TeV LHC is easily calcu- 
lated 0, and this falls from a few pb at mz ~ 300 GeV 



to about 10 ab for 



1.5 TeV. As a typical estimate 



of the integrated luminosity of the LHC is about 3 ab 



the latter value corresponds to around 30 fifi pairs — 
which may be just enough for a discovery. Accordingly, 
we take m^^ < 1.5 TeV as a reasonable criterion for a 
possible LHC discovery of the light stop. 

This light stop criterion is not, however, the only re- 
striction on the parameter space of the SUSY model in 
question. The requirement of a light Higgs state in the 
range 122 — 129 GeV, coupled with the recent observation 
of the rare decay Bg — )■ M~^J4~ @i the updated B ^ tv 
branching ratio at BELLE Q , and the negative results of 
all direct searches for new particles at the LHC [1] , taken 
together, make a formidable combination, imposing tight 
restrictions on models of new physics. Among the most 
severely affected are models where the electrowcak sector 
differs from that of the SM. All SUSY models fall in this 
class. In an article [1] written last summer, just before 
the announcement of the new boson discovery, two of the 
present authors had shown that even the most restrictive 
model of supersymmetry, viz. the constrained Minimal 
Supersymmetric SM (cMSSM) was affected by the 
then-available data only marginally — at the periphery 
of its parameter space — leaving most of its parameter 
choices viable. The new data obtained since then have, 
in the span of a few months, changed the situation quite 
dramatically — the allowed parameter space now appears 
to be severely restricted [lOt . 

Since the new data are all compatible with the SM, 
they may be expected to drive SUSY towards the de- 
coupling limit. This, of course, runs contrary to the 
light stop scenario desired by us. In the current arti- 
cle, we examine whether the cMSSM can accommodate 
these conflicting requirements. Moreover, we also con- 
sider a popular extension of the cMSSM, known as the 
Non-Universal Higgs Mass (NUHM) model [ll|, which 
could, in principle, evade some of the constraints which 
affect the cMSSM adversely. Both the models in question 
have been widely discussed in the literature [H, [ill ' s-'^d 
do not require to be extensively reviewed in a focussed 
study such as the present one. Below, we just summarize 
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the principal inotivation(s) and the parameter choices. 

The cMSSM is a theoretically well-motivated super- 
symmetric model with the least number of parameters, 
and consequently has the highest predictivity of all low- 
energy SUSY models. The free parameters of this model 
are (i) the universal fermion mass mij^ at the grand uni- 
fication (GUT) scale, (ii) the universal trilinear coupling 
at the GUT scale, (iii) the ratio tan /3 of vacuum ex- 
pectation values of the two Higgs doublets at the elcc- 
troweak scale, and (iv) the universal scalar mass mg at 
the GUT scale. In addition, we can choose the sign of 
the Higgs mixing parameter /i to be positive or negative. 

Predictivity is the most appealing feature of the 
cMSSM, though it also makes it more vulnerable to all 
sorts of experimental constraints. Therefore, even if the 
cMSSM is ruled out by the data — or, more likely, pushed 
to its decoupling limit — it would be logical to look for its 
close cousins, i.e. models which contain at least some of 
the more attractive features of the cMSSM. The NUHM 
model is perhaps the most natural of these choices, for 
all that it does is to relax one of the cMSSM conditions 
and make the two Higgs mass parameters M^.^ and M^,^ 
different from the universal scalar mass toq at the GUT 
scale. 

In our numerical studies, we study the independent 
variation of these parameters in the ranges 



TOi/2 6 [0, 2] TeV, 
Ao e [-10, 10] TeV, 
Mk, e [0, 5] TeV, 



mo G [0, 5] TeV, 
tan;3 G [5, 60], 
M,,, G [0, 5] TeV, 



which are compatible with theoretical considerations and 
cover the entire region which allows for a light stop <i of 
mass below 1.5 TeV. We also restrict ourselves to ^ > 0, 
since only the positive sign of can ameliorate the 3.6(j 
disagreement of the SM with the measurement of the 
anomalous magnetic moment {g — 2)^ of the muon Q. 
However, as it is still not possible to be compatible with 
the {g — 2)^ measurement to better than 95% C.L. [l^ . 
we choose to be conservative and do not impose this par- 
ticular constraint in any other way. Having chosen the 
parameter ranges, then, we make a random scan over 
the entire range of these variables. For each set of input 
parameters, we use the public domain software SUSY- 
HIT [15| to compute the values of the entire supersym- 



metric spectrum and couplings at the respective scales 
at which the experimental data are available. These are 
then taken as inputs to calculate a slew of high- and low- 
energy observables for which experimental data are avail- 
able. For this part, we use the software SuperISO [T6j 
and FeynHiggs [rTj . 

In order to constrain the parameter space, we first im- 
pose the "theory" constraints, i.e. the obvious require- 
ment that the electroweak symmetry be broken at pre- 
cisely the electroweak scale, and the requirement that the 



lightest sparticle be neutral (which arises from its iden- 
tification dark matter candidate). However, tak- 
ing a conservative stance again, we do not impose con- 
straints arising from the relic density of this neutralino, 
since such constraints can be evaded in scenarios with 
non-thermal dark matter [l8| . In addition, we use the 
following constraints, many of which have been recently 
updated (mostly during 2012): 

(i) The mass of the lightest Higgs boson Mh lies in the 
narrow range 122.5 GeV < Mh < 129.5 GeV. This range 
includes the 95% C.L. errors from the ATLAS and CMS 
data [l| as well as a 1.5 GeV theoretical uncertainty [l9j . 

(ii) The branching ratios B for some rare B decays are 
constrained as follows, at 95 % C.L. llllUllI 



1.1 X 10~^ < B{Bs^ fi-^fi') < 6.4 X 10"^ , 
B [Ba fi+fi') < 9.4 X 10"^° , 
2.6 X lO"'' < B{B ^ Xs-/) < 4.5 x lO"'^ , 
0.44 X 10""^ < B {B TVr) < 1-48 x lO"'' . 

As a matter of fact, the first measurement of B{Bs /^/i) 
has recently been reported at the LHCb Also, the 
measurement of B{B — >■ tv) with hadronic tag has re- 
cently been updated by BELLE Q and brings the value 
of this ratio to a much smaller value than earlier mea- 
sured. This not only relaxes the tension of this mea- 
surement with the SM prediction, but also dissipates the 
strong constraints it had put on the cMSS M [22| . Here 
we take the BELLE average of B{B tv) j2lj . 

(iii) Another low energy measurement arises from the 
2006 measurement, at the DA$NE facility at Frascati, 
of the process — ?> jJi^v^. Taking the 95% C.L. limits, 
this gives us [1^ 

0.6331 < B{K+ fi+i^^) < 0.6401 , 

which can be very restrictive for a certain class of elec- 
troweak models. 

(iv) We also impose the 95% C.L. lower bounds on the 
masses of sparticles, as listed by the Particle Data Group 
(PDG) Q. These are given below, with all masses being 
measured in GeV: 



Mc;0 > 46 

Al 



rUg > 1100 



mj> 107 . 



M^o > 116 , mf^ > 95.7 , m,7f > 94 , 
M~± > 94 , m^^ > 89 , to?^ > 81.9 
Mg > 500 . 



Here q and vi refer to superpartncrs of the first two 
generations. There are more specific constraints on the 
gluino mass Mg, which depend on the mass m^^ of the 
light stop [1] ; we have verified that these are satisfied for 
our allowed points. 

Note that in selecting the allowed points in our ran- 
dom scan, we do not treat the early measurements of the 
Higgs boson branching ratios as constraints, since these 
are based on a rather small sample of Higgs boson decay 
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FIG. 1: Allowed regions in the parameter space for fi > 0. We consider two sections of the parameter space, viz. (a) the 
mi/2-mo plane, and the (b) y4o-tan/3 plane. The black dots indicate the cMSSM, while the green (grey) ones indicate the 
NUHM. 



events and, consequently, have large error bars. Instead, 
we shall presently argue that these measurements could 
constitute stringent tests for the models in question as 
more data become available. 

Taking all these constraints into account, the allowed 
points in our random scan are shown in Figure [1] The 
left panel — marked (a) — shows the mo-mi/2 plane, 
while the panel on the right — marked (b) — shows the 
Ao-tan/3 plane. We should note that all the parame- 
ters vary in these plots, 6.17. in the panel (a), ^0 a-nd 
tan /? arc not held fixed but are allowed to vary over the 
full ranges mentioned above. A glance at the panel (a) 
shows that the two models, viz. cMSSM and NUHM, 
hardly differ so far as the allowed region is concerned. 
Only in the so-called stau co-annihilation region, where 
mo is small and mi/2 is large, does the NUHM enjoy 
a slightly more viable status than the cMSSM. However 
— and this is mostly due to the Higgs boson mass mea- 
surement — both the models are severely constrained in 
the low- mo and low-mi/2 region. Unfortunately, this is 
also the region in which large signals for supcrpartners 
are predicted at collider machines, such as the LHC, so 
this result is unfavourable to the school of thought that 
supersymmetry is "around the corner" . The paucity of 
points in the "decoupling" region of large-mo and large- 
™i/2 is a direct consequence of imposing the restriction 
"^ti < 1-5 TcV. However, the large number of allowed 
points in the central region of the panel clearly indi- 
cates that these two models — with all their economy 
of parameter choices — are still capable of explaining all 
the experimental data, as well as predicting a light stop. 
This conclusion is also reflected in the panel (b), where 
it is clear that large negative values of Aq are favored, 
but otherwise the only real restrictions are on large val- 
ues of tan/3. These arise from the tan^ /3 sensitivity of 



B{Bs 01, coupled with the regime imposed 

by our artificial cutoff m^^ < 1.5 TeV. In the cMSSM, 
this forces us to have tan/3;$45, but the NUHM can still 
have allowed points with larger tan /3 values. 

It is interesting to ask why large negative values of 
Aq are preferred. This is because the measured mass 
of the Higgs boson Alh — 125-126 GeV is sufficiently 
removed from the decoupling limit — typically quoted Q 
as Mh w 119 GeV for ^0 = — to require a light stop 
ti to explain that deviation. In order to sustain a light 
stop, while simultaneously keeping the squarks of the first 
two generations heavy, we require a large "seesaw" effect 
in the stop mixing matrix, and this happens when large 
negative values of Aq appear in the off-diagonal terms of 
that matrix. 

This is nicely illustrated in Fig. [21 where we have shown 
the scatter of parameter choices with the mass m^^ plot- 
ted against the Aq parameter, and the effect of different 
classes of constraints on the allowed region in this plot. 
In panel (a), we impose only "theory" constraints, which 
are not only indifferent to the sign of Aq, but also al- 
low a stop lighter than 100 GeV in the Ao « region. 
In panel (b), we impose, in addition, the indirect con- 
straints arising from low-energy measurements, includ- 
ing B-decays, and this immediately affects the Ao ~ 
region, driving the stop mass above 120 GeV. Imposition 
of the direct search constraints, especially the require- 
ment that rriq > 1.1 TeV for the first two generations [1], 
cleans up the Aq ~ region in panel (c) much more ef- 
ficiently, though a few points with a stop mass between 
120-200 GeV still survive. Here, too, both signs of Ao are 
equally preferred. However, when we impose the Higgs 
boson mass constraint, the allowed parameter space in 
panel (d) is affected quite severely. Not only arc small 
values of ^0 completely washed-out, with Aq being driven 
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FIG. 2: Illustrating the correlation between the mass of the 
light stop and the Aq parameter, as different sets of con- 
straints are imposed successively. Each panel includes the 
constraints from the previous panel. The shaded region repre- 
sents the LEP constraint on mi . Other conventions followed 



are exactly the same as in Fig. |l] 



below —2 TeV or above 8 TeV, but the stop mass is es- 
sentially pushed above 500 GeV, with only a few outlying 
points in the range 250-500 GeV (which may very well 
be due to specially fine-tuned cancellations). Thus there 
is no real surprise in the fact that experimental searches 
for stops lighter than 550 GeV at the 7-8 TeV LHC have 
come up empty-handed [2^. On the other hand, it is 
encouraging that there is a profusion of allowed points 
with a ti mass in the range 600 GeV-1 TeV — which, 
as we have argued, could easily turn out to be a happy 
hunting ground at the 14 TcV LHC. 

We now come to the final test, viz, whether the mod- 
els in question can predict the observed branching ratios 
of the Hi ggs boson as observed in the ATLAS l2g and 
CMS |27H3G | data. Of course, what is observed at the 
LHC is not the branching ratio per se, but a number of 
events of a certain kind which can be identified as arising 
from the decay of a Higgs boson resonance. This cor- 
responds to the convolution of the integrated luminosity 
C with a production cross-section and a branching ratio, 
i.e. for a particle pair XX, we have 



We then define the so-called signal strength 



XX) 
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(ysuipp^hO) Bsm(/i°^AA) 



(1) 



, (2) 



where the label SM indicates the SM predictions with 
nih = 126 GeV. The Higgs production cross-sections are 
calculated using FeynHiggs and the branching ratios 
are calculated using the software HDecay, which is in- 
corporated in the SUSYHIT package and is used in con- 
junction with the supersymmetry parameters allowed by 
the constraints described above. Our results can then be 
compared with the latest experimental values of U y y 
given by the two experimental collaborations |26l - l30j . 
which are listed in Table |T] below. 



ATLAS 



CMS 



fJ-WW 


L0±0.3 \26} 


0.76 ±0.21 [27| 


fJ-zz 


1.5 ±0.4 [26| 


o.9ii;;-i° [28] 




0.8 ±0.7 [26| 


1.1 ±0.4 [29| 




1.6 ±0.3 \26} 


0.78+0:^6 [30] 



TABLE I: Experimental data on Higgs boson signal strengths, 
as reported in March 2013 by the ATLAS and CMS [13- 
[30I ] Collaborations. For fj,-yj we use the multivariate analysis 
result from CMS 

From Table U it is clear that at 95% C.L. all the mea- 
surements are compatible with a SM Higgs boson. We 
also note that these results are based on almost all the 
data collected in the 7-8 TeV run of the LHC, and hence, 
no substantial improvement may be expected till the en- 
ergy/liminosity upgrade is in place. 
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FIG. 3: Correlation plots for the signal strengths of the Higgs boson in different channels, vis-a-vis ATLAS and CMS bounds. 
SUSY conventions follow those of Fig[T] In each panel, the SM prediction lies at the crossing of the two short (orange) lines, 
which denote the theoretical uncertainty due to the spread in Higgs boson mass. 



With these circumstances in mind, we now make a 
comparison between the signal strengths predicted in the 
cMSSM and the NUHM versus the experimental data. 
These results are shown in Fig. [3] In the panel on the left 
— marked (a) — we plot fJ-ww versus ^zz, with black 
dots indicating the cMSSM and green (grey) dots indi- 
cating the NUHM. The blue (solid line) boxes indicate 
the ATLAS limits at 1(t and 2a respectively, while the 
red (dashed line) boxes indicate the CMS limits, again at 
Icr and 2(7 respectively. The SM prediction lies at (1.0, 
1.0), i.e. at the crossing of the two short (orange/grey) 
lines, whose lengths correspond to the uncertainties aris- 
ing from the allowed range in Higgs boson mass. It is im- 
mediately obvious that the CMS data are slightly more 
restrictive than the ATLAS data. If we take the bounds 
at 2(7, then it is obvious that the full scatter of points in 
both the cMSSM and NUHM are allowed by the data, as 
is, of course, the SM. On the other hand, if we consider 
the 1(7 constraints, then the SM is disfavoured, whereas 
the two supcrsymmetric models survive. 

The comparison changes a little if we turn to the panel 
on the right — marked (b) . The conventions for this plot 
are exactly the same as before, the only difference being 
that we have plotted the more controversial versus 
firr, the former being highly sensitive to new physics con- 
tributions. Here the ATLAS and CMS results are some- 
what different. The SM point as well as the scatter of 
points in the two supcrsymmetric models lie in the heart 
of the CMS-allowed region, whereas for the ATLAS re- 
sult, the SM lies at the 2a border, while most of the su- 
pcrsymmetric points lie further out. Thus, in the ATLAS 
data, at least, there is some hint of new physics, but the 
CMS data appears to confound any such speculations. 

It is quite clear from this study of the signal strengths 
that — as of now — both the cMSSM and NUHM survive 



all the experimental tests, even if we demand a light stop 
state discoverable at the LHC. In fact, they do as well, 
if not better, than the SM — as is only to be expected 
from models with more free parameters. The predictions 
of these models in their allowed ranges may deviate sig- 
nificantly from the SM predictions; however the small- 
ness of the scatter compared to the experimental errors 
in both the plots in Fig. |3] shows that, the present status 
of these models is not so different from that of the SM. 
Indeed; in both the panels of Fig. |3] the SM point lies 
right in the middle of a dense clustering of allowed SUSY 
points. This indicates that these two models can mimic 
the SM predictions for these signal strengths even with 
a light stop of mass m^^ < 1.5 TeV. It is thus much too 
early to write off any model of supersymmetry, even the 
highly restrictive cMSSM. 

As regards the future, one can only speculate. Obvi- 
ously, when the results from the analysis of more data 
become available, the boxes in Fig. [3] will shrink, corre- 
sponding to a reduction in the experimental errors. Such 
changes are often accompanied by changes in the cen- 
tral value(s) as well, which would appear on the plots 
as a translation of the box(es) in some direction. Three 
distinct possibilities may now be considered. 

(i) The new boxes may zoom in to the SM prediction. 
In this case our plots in Fig. [3] show that there are always 
points in the parameter space of the cMSSM and NUHM 
which will give the same fixx predictions as the SM. 
Thus, these models will survive — and even predict a 
discoverable light stop — though with even more severe 
constraints on the parameter space. 

(ii) The shrinking boxes may move away from the SM 
predictions, but stay within the scatter of allowed SUSY 
points shown in Fig. |31 This, would, of course, form a 
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powerful argument for these minimal versions of SUSY 
and would certainly enhance our expectations for a light 
stop discovery. 

(iii) In addition, we must also consider the possibility 
that the new data will disfavour both the SM and these 
two highly-predictive SUSY models — in which case one 
would have to invoke less restrictive versions of SUSY 
(e.g. the 19-parameter pMSSM [T^]) or other forms of 
new physics to explain the results. 

Before concluding, we briefly consider the case when 
the "under the lamp post" restriction of m^^ < 1.5 TcV 
is removed, or rendered irrelevant by negative results of 
direct stop searches at some time in the future. We 
find that relaxation of this requirement does not affect 
the low- mo and 1ow-TOi/2 region, but does populate the 
high-mo and high-7711/2 "decoupling" region. At the same 
time, high values of tan /3 are marginally less constrained. 
Moreover, small negative values of Aq are now allowed, 
but we should note that these correspond only to heavier 
values of m^^ > 1.5 TeV. More interestingly, the scatter 
of points in both panels of Fig. [3] does not change ap- 
preciably. It follows that our remarks about the future 
of tests involving Higgs boson decays would generally re- 



main valid even with a heavy m^^ , which may be beyond 
the kinematic reach of the LHC. 

To sum up, then, we have made a critical study of 
two of the most restrictive models of supersymmetry, 
cMSSM and NUHM, taking into account the recent mea- 
surements of the Higgs mass and the updates in many 
flavour-physics observables. We find that these models 
can actually survive all the constraints and predict a light 
stop which could be discovered, in principle, at the LHC. 
The allowed parameter ranges for both the models seem 
almost identical, except for tan/3, which is restricted to 
be ^45 in the cMSSM with a light stop, while it could be 
larger in the NUHM. The branching ratios of the Higgs 
boson to WW, ZZ, tt or 77 can be an important testing 
ground for the future of these models. However, as long 
as the experimental data do not exhibit definite disagree- 
ment with the SM, these SUSY models will continue to 
survive, and may even offer the exciting possibility of a 
light stop discovery in the coming runs of the LHC. 
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